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We study the singlet part of the free energy of a static quark anti-quark (QQ) pair 
at finite temperature in three flavor QCD with degenerate quark masses using Nt = 
4 and 6 lattices with Asqtad staggered fermion action. We look at thermodynamics 
of the system around phase transition and study its scaling with lattice spacing 
and quark masses. 



^ ■ 1. Simulations and Results 

Static free energy of a quark- antiquark pair is a frequently used tool to 
1^ ■ study non-perturbatively the in-medium modification of inter-quark forces. 

' It is calculated as the difference in the free energy of the system with static 

quarks and the same system without them, while the temperature remains 
. constant. Due to colour symmetry such quantity contains the singlet and 

^ I the octet contributions, and the usual definition of the free energy is thus 

referred to as the colour-averaged. This quantity has been extensively in 
SU{2) and S'[/(3) gauge theories e.g. [3] and full QCD [8]. However, singlet 
and octet channels were considered in detail only in pure gauge theory [2, 5, 
6]. In the case of full QCD the first results for singlet and octet free energy 
for two flavor QCD have appeared only very recently [9]. Here we present 
our results for 3 flavor QCD using the so-called Asqtad staggered fermion 
action [10] with two different lattice spacings (corresponding to A^t = 4 and 
6) at three different quark masses. 

Our analysis is to a large extent based on the gauge configurations 
generated by the MILC collaboration using the Asqtad action. Therefore 
we adopt their strategy for fixing the parameters which is described in 
Ref. 11. We use the most recent value of ri extrapolated to continuum and 
to the physical value of the light quark masses ri — 0.317 fm to convert the 
lattice units to temperature. 

The free energy of a static quark anti-quark pair contains a lattice spac- 
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ing dependent divergent piece, and thus needs to be renormalized. Follow- 
ing Ref. 5 we do so by normalizing it to the zero temperature potential at 
short distances where the temperature dependence of the free energy can be 
neglected. The static quark potential has been studied by the MILC collab- 
oration at three different lattice spacings and various quark masses [11]. Wc 
use the following form of the zero temperature potential, which reproduces 
MILC data well in the whole range of masses and lattice spacings 
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In Fig. 1 we show that both the potential and effective coupling constant 
as{r) defined as as(r) = 3/4r^rfF(r)/rfr are reproduced well. Results pre- 
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Figure 1. Zero-temperature static potential (left) and the running coupling (right) [data 
points from MILC collaboration] 



sented here come from MILC lattices as well as our simulation for quark 
masses mught = 0.2ms, 0.4ms, 0.6ms on lattices 12^^ x 4, 8"^ x 4 and 12"^ x 6. 
The temperature range was 135 — 412 MeV for the first two lattice sizes 
and 145 - 310 MeV for the last one. 

Following Ref. 14 the free energy of static quark- antiquark ( QQ) pair 
in the color singlet channels is defined as 

exp(-Fi (r, T)/T + C) = ^Tr(W{r)W^ (0)) (2) 

Here W{x) = H^^o ^ Uq{t,x) is the temporal Wilson line, L{x) = TrW{x) 
is known as the Polyakov loop. As W{x) is not gauge invariant one needs 
to fix a gauge. Here we will use the Coulomb gauge as advocated in [4]. 
This approach is exactly valid at zero temperature and is numerically true 
at finite temperature. One can also consider the color averaged free energy 
defined as 

oM-Fav{r,T)/T+C) = i(L(f)Lt(0)) (3) 
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We start the discussion of our numerical results with the case of the 
quark of mass 0.4ms on the 12'^ x 4 lattice. The corresponding numerical 
results for the singlet free energy are shown in Fig. 2. The free energy 
approaches a finite value Fl^{T) = limr^oo FiifjT), i = l,av at large 
distances which is visually interpreted as string breaking at low temperature 
and screening at high ones. Note that the distance where the free energy 
effectively flattens is temperature dependent, it becomes smaller at higher 
temperatures. At small distances the singlet free energy is temperature 
independent and coincides with the zero temperature potential, as at small 
distances medium effects are not important. The numerical results for other 
values of the quark masses are similar. The scaling with the lattice spacing 
is remarkably good. 




Figure 2. Singlet free energies and the effective screening radius 



To characterize the range of interaction in the medium it is convenient 
to introduce the effective screening radius Vgcr defined as: Fi(r = r^cr, T) = 
0.9F^{T). Here F};^{T) is the asymptotic value of the singlet free energy 
at infinite separation. In Fig. 2. we show the values of Vscr for three differ- 
ent quark masses and 12^ x 4 lattices. Certainly as Fi(r, T) has statistical 
errors it is difficult to determine at exactly which distance r the equation 
Fi{r = TscriF) = 0.9F^{T) holds. We have tried to estimate this un- 
certainty in the values of rscr and show them in Fig. 2 as errors bars. At 
small temperatures the value of the screening radius is about 0.9/m and is 
temperature independent. As we increase the temperature Vscr decreases 
reaching the value of 0.5/m at the highest temperature. Note that the 
temperature dependence of rscr is roughly the same for all quark masses. 
On Fig. 3 (left) we plot the asymptotic value of the free energies; the quark 
mass dependence likely vanishes at small temperatures (T < 150MeV) and 
definitely negligible at high temperatures (T > 250MeV). It is however 
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Figure 3. Infinite separation free energies for various quark masses (left) and internal 
energy at 0.4ms 



significant in the transition region. 

In the hmit of very small temperatures we expect Foe (T) to be temper- 
ature independent and related to twice the binding energy of a heavy-light 
{D— or -B— ) meson 
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Based on this observation in Ref. 15 it has been argued that the decrease 
of Fca{T) with the temperature close to Tc implies the decrease of the 
Mo.B leading to quarkonium suppression. However, Foo{T) also contains 
an entropy contribution due to the presence of a static QQ pair: 

dF^{T) 



SooiT) 



dT 



(5) 



Then we can calculate the energy induced by a static quark-anti-quark pair 



(6) 



Numerically the derivative with respect to the temperature in Eq. (5) was 
estimated using forward differences. 

On the right side of Fig. 3 we show the energy 1/^0 as function of tem- 
perature. Both the entropy and the energy show a strong increase near 
Tc- This large increase in entropy and energy is probably due to many- 
body effects and makes the interpretation of Uoo as the binding energy of 
heavy- light meson not very plausible. 



2. Conclusions 

The free energy gets screened beyond some distance for all temperatures 
as expected. For small temperature this distance, the effective screening 
radius, does not depend on the temperature and is about 0.9/m. As the 
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temperature increases the effective screening radius decreases. Light quark 
mass dependence of the screening radius is neghgible within our statistical 
accuracy. We have also identified the entropy contribution to the free energy 
as well as the internal energy at large distances and found that they show 
strong increase at T^. We have found substantial quark mass dependence 
of the free energy in the vicinity of the transition. 
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